Neuronal networks with coupling through amyloid beta: towards a theory for 

Alzheimer's disease. 
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Alzheimer's disease (AD) is a common form of dementia observed in the elderly due to neurode- 
generative disorder and dysfunction. This arises from alterations in synaptic functioning of neurons 
leading to cognitive impairment and memory loss. Recent experimental studies indicate that the 
amyloid beta (A/3) protein, in its dimer and oligomer forms, affects the synaptic activity of neurons 
in the early stages of AD. However, the precise mechanism underlying A/3 induced synaptic depres- 
sion is still not clearly understood. In this paper, we introduce an electrical model that provides a 
possible mechanism for this. Our studies show that the competing effects of synaptic activity and 
the indirect interaction mediated by A/3 , can disrupt the synchrony among neurons and severely 
affect the neuronal activity. This then leads to sub-threshold activity or synaptic silencing. This 
is in agreement with the reported disruption of cortical synaptic integration in the presence of A/3 
in transgenic mice. We suggest that direct electrophysiological measurements of A/3 activity can 
establish its role in AD and the possible revival, proposed in our model. The mechanism proposed 
here is quite general and could account for the role of the relevant protein in other neuronal disorders 
also. 
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INTRODUCTION 



Alzheimer's disease (AD) is the most common cause 
for dementia characterized by a progressive neuronal dys- 
function leading to neuronal death and depressed brain 
function [l|-[3|. It has been established over the past 
decade that pathologically this neurodegenerative disor- 
der may arise from extra cellular deposits of Amyloid /3 
(A/3) protein in senile plaques and inter neuronal neuro- 
fibrillatory tangles However, recent experimental 

studies indicate that the amyloid beta (A/3) protein, in 
its dimer and oligomer forms, affects the synaptic ac- 
tivity of neurons in the early stages of AD [5-8]. The 
actual mechanism by which A/3 oligomers participates in 
the impairment of synaptic function is still not clearly 
established. 

In this work, we introduce a physical model that can 
explain the mechanism of A/3 induced synaptic depres- 
sion especially in the early stages and point towards pos- 
sible methods for prevention or reversal. We start with 
the simple hypothesis that since the primary functions of 
neurons are through electrical activity, it is quite possible 
that the disturbance or impairment of this function is also 
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electrical in nature. This is supported by the fact there is 
no evidence of physical blocking or damage of synapses or 
neurons in the very early stages [7] . Amyloid beta, being 
a protein is not expected to have any electrical activity of 
its own and no such activity has been reported so far. But 
there exists a possibility of an induced electrical activity 
(for example, through polarization) caused by electrical 
signals generated by nearby neurons. Recent molecular 
dynamics based study on A/3 indicates the possibility of 
such an induced dipole moment due to rearrangement of 
charges in the presence of an electric field [9[. This in- 
duced electrical activity of A/3, in turn, provides an input 
back to the surrounding neurons (giving rise to a feed- 
back loop). Thus A/3 can act as an intermediate agency 
that provides an extra channel of interaction among neu- 
rons. This is in addition to the regular synaptic interac- 
tion between the neurons. The normal coupling through 
synapse leads to in-phase synchronization of the activ- 
ities of neurons. The additional coupling through A/3 
is equivalent to an environmental coupling of the type 
discussed in a previous paper where it was shown 
that it can lead to an anti-phase synchronization. Thus 
the above two types of couplings between the neurons 
represent two competing synchronizing tendencies. Both 
these tendencies become effective above their respective 
thresholds for the coupling strengths. In addition, if their 
effects balance each other, we expect the electrical activ- 
ity of the neurons to be severely affected. The numerical 
simulation carried out on a neuronal network with ad- 
ditional coupling through A/3 establishes disruption of 
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activity and silencing. 



2. NEURONS COUPLED INDIRECTLY VIA A/3 

The mechanism proposed in the paper is appHed to a 
network of 10 neurons using the Hindmarsh-Rose ( HR) 
model [10]. A simple way of modeling the effect of the 
extra channel via A/3, in such a network is given below. 
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(1) 



The first three equations in ([T]) come from the HR model 
(except for the w term) and represent the neuronal ac- 
tivity in a normal brain. Here the variable x represents 
the voltage of the neuron and y and z are variables re- 
lated to ion currents and the neurons are coupled via ex- 
citatory synaptic couplings [21] represented by the last 
term in the first equation. We model the pathologi- 
cal activity of A/3 in the case of AD by the additional 
variable w. The nature of the electrical activity of A/3 
considered here is not inherent but induced. So with- 
out input from neurons, this induced electrical activity 
should decay. This is represented by the decay term, - 
KW^ in the w equation where hz is the damping parame- 
ter. The second term in the w equation represents the 
input from the surrounding neurons. The additional CaW 
term in x equation is the extra input to the neuron volt- 
age through A/3 as induced electrical activity. Here, we 
choose the parameters a = 3,6 = 5, p = 0.006, 5 = 4, 
X = 1.6,/ = 3.2, such that the individual neurons have 
chaotic responses. The damping parameter of the amy- 
loid beta is taken to be z^: = 1 . The synaptic coupling is 
modeled by the sigmoidal function given in the last term 
in the first equation. Here the reversal potential Vr < Xi 
makes the synapse excitatory. The other parameters in 
the synaptic coupling term are chosen to be ^ = —0.25 
and A = 10. The network is constructed using an adja- 
cency matrix A, whose elements are randomly chosen as 
either or 1 in Eq.(l). Aij = 1 indicates that i^^ neuron 
receives synaptic input from j^^ neuron. Here the matrix 
A is asymmetric, such that the coupling is mostly unidi- 
rectional though bidirectional couplings can also occur. 
The synaptic input represented by a sigmoid function is 
the self-feedback current of i^^ neuron modulated by the 
j^^ neuron. Using random initial condtions, the network 
is evolved using Runge-Kutta algorithm with step size 
0.01 for sufficiently long time. 

First we study the baseline activity by considering a 
network of 10 neurons with HR model as the nodal dy- 
namics and coupled via excitatory synapses alone ( ig- 
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FIG. 1: Time series of the voltage of all neurons (grey) and 
their average (black). Here a random network of 10 neurons 
is used with different strengths of synaptic coupling in the 
absence of coupling through A/3. (A) Unsynchronized state 
for Cs = 0.01 . Here the average voltage is much smaller than 
individual spikes. (B) Synchronized bursting for Cs = 0.5 and 
(C), Synchronized states for Cs = 1.0. HR parameters chosen 
are as a = 3,6 — 5,p — 0.006,s — 4,x — 1.6 and / — 3.2, for 
which individual neurons are in the chaotic state. 



noring the influence of A/3). For very small strengths of 
synaptic coupling , the neurons do not act in synchrony. 
This is clear from the sequences of neuronal voltages plot- 
ted in Fig. [TJA. Here the grey lines give the neuronal 
voltages from all the neurons while the black shows the 
average response. The emergence of synchronous behav- 
ior becomes evident as the synaptic strength increases. 
Thus, Figs.[TjB and[TjC correspond to two typical cases 
of coordinated activity when only synaptic interaction is 
present. 

Next, we introduce the extra coupling through A/3 in 
the same network and the consequent changes in neu- 
ronal voltages are shown in Fig. 2.A-F. We see that the 
competing effects of the two modes of interaction lead to 
severe disruptions of the neuronal activity which result 
in sub threshold behavior and in some cases to complete 
silencing. Further, the synaptic activity at one neuron is 
calculated using the last term in Eq. ([1]) and is shown in 
Fig. [2lG,H,I. In Fig. [2l we plot the total synaptic activ- 
ity at one neuron and for different values of A/3 coupling 
,show that A/3 coupling also disrupts the synaptic behav- 
ior (Fig.2.G-I). We note that such disruptions in synaptic 
behavior due to A/3 have been reported in earlier exper- 
imental observations [l2|. The induced activity of A/3 
based on our mechanism is evident from the behavior of 
the variable w shown in Fig 2.J-L. This induced activity 
of A/3 is caused by the electrical signal from the neurons 
and when it is of sufficient strength, it can cause neuronal 
dysfunction. We also observe that the route to synaptic 
silencing or sub threshold activity is through a stage of 
desynchronization of the neuronal outputs. This is clear 
from Fig.2.B and E where the strength of is just not 
sufficient for synaptic depression. We note that this is 
consistent with the measurements of cortical neuronal 
signals reported in transgenic mice [l3|. 

That the resultant behavior of neurons under the 
two competing interactions depends on their relative 
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FIG. 2: Disruption of neuronal activity by A/3. (A)-(F) Time series of the voltage of all neurons in the network (grey) 
and their average (black) for different strengths of synaptic coupling and coupling through A/3. (A) Normal synchronized state 
(cs = 0.50, ea = 0.01) (B) Desynchronized state {cs = 0.50, Ca = 0.50). (C) Sub-threshold activity (e^ = 0.50, Ca = 3.50) 
(D) Normal synchronized state (e^ = 1.00, ea = 0.01). (E) Desynchronized state (cs = 1.0, ea = 1.60) (F) synaptic silencing 
(es = 1.00, Ca = 2.00). (G-I) Synaptic input to one neuron and (J-L) Induced activity on A/3 as a function of time at different 
strengths of synaptic coupling and coupling through A/3 (G) Synaptic input for synchronized state (e^ = 1.00, Ca = 0.01) (H) 
synaptic input for desynchronized activity near sub-threshold oscillatory state ( es = 0.50, ea = 0.50 ). (I) synaptic input for 
desynchronized activity near neuronal silencing (es = 1.0, Ca = 1.60) (J) Induced activity on A/3 due to signal from neurons for 
Cs — 1.00, Ca — 0.01 (K) Induced activity on A/3 for ts — 0.50, Ca — 0.50. (L) Induced activity on A/3 for ts — 1.0, Ca — 1.60. 



strengths is evident from Fig. 2. Hence we study in de- 
tail the neuronal activity in the parameter plane of the 
coupling strengths (ea, ). The maximum value of the 
voltage of neurons in the network is used as an indicator 
to identify regions of synaptic silencing in the parame- 
ter plane of coupling strengths. This is shown in Fig. [3l 
where black regions correspond to complete electrical in- 
activity (quiescent state) and red regions correspond to 
sub-threshold oscillations. We find two patches of red 
and black regions in the parameter plane where the two 
competing effects balance each other. 

Based on our model, we propose that one of the possi- 
ble mechanisms to prevent this regulatory activity of A/3 
is either inhibiting electrical activity on A/3 or by reduc- 
ing the value of parameter n which amounts to very slow 
relaxation of the induced activity of A/3. The latter case 
is shown in Fig. 4, where by decreasing the synaptic 
activity is found to be restored. 



3. DISCUSSION 

We have introduced a physical mechanism that ex- 
plains the electrophysiological death of neurons due to 
the influence of A/3 oligomers in the context of AD. We 
note that some of the key features of our analysis agree 
with the findings from other related and direct obser- 
vational studies. Thus it is found by electrophysiological 
measurements in cultured hippocampal neurons that ele- 
vated A/3 production disrupts both presynaptic and post- 
synaptic functions Our results also indicate that 



synaptic depression occurs in all the neurons. A study of 
hippocampal neurons reported that the synaptic depres- 
sion requires neuronal activity [l5|. In our model also, 
with no neuronal activity, A/3 will not become electrically 
active and therefore cannot cause depression. We observe 
that the synaptic silencing is associated with considerable 
disruption of synchronized activity. This is in agreement 
with the reported disruption of cortical synaptic integra- 
tion in the presence of A/3 in transgenic mice [13]. 

Some of the mechanisms for AD reported recently 
are either biochemical due to generation of free radicals 



16 1 or chemical due to binding of A/3 in cell membrane 
17,[l8|. The mechanism introduced in this paper is novel 



in that it is the first electrical model to the best of our 
knowledge. Moreover, our model for the first time pro- 
poses the presence of induced activity in A/3 oligomers 
when synaptic disruption occurs. We note that indirect 
evidence for induced electrical activity can be drawn from 
the reported studies on birefringence in certain amyloid 
proteins such as /3-Lactoglobulin under applied electrical 
pulses [19]. This indicates the possibility of such an ac- 
tivity in A/3 under the influence of neuronal signals. A 
direct test for such induced activity can be obtained by 
electrophysiological measurements in cultured media or 
brain slices. Such experiments can also lead to possible 
revival strategies by external interventions of the type 
proposed in our model. 

None of the studies reported so far explain precisely 
the reason why A/3 monomer is not effective in inducing 
synaptic suppression but dimers and other n-oligomers 
are. We conjecture that a possible explanation through 
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our model would be that the monomer does not develop 
electrical activity of the required magnitude to have a 
measurable effect, while higher oligomers develop electri- 
cal activity of the same order of magnitude as a normal 
synaptic interaction. A direct measurement of the in- 
duced effect on amyloid monomer and dimer in a field 
equivalent to that across the synapse can give direct sup- 
port for this conjecture. Thus the theory proposed here 
opens up interesting observational research in at least 
two directions. We have obtained qualitatively similar 
results with other models such as Hodgkin-Huxley and 
FitzHugh- Nagumo models. These results are presented 
in the supporting information. Finally, we propose that 
the mechanism introduced here is quite general and can 
be applied in the context of other neurodegenerative dis- 
eases such as Parkinson's disease, Huntington's disease, 
Creutzfeld- Jakob disease, etc where specific protein de- 
posits have been identified as the neurotoxic agents [2Q|. 



FIG. 3: Regions of disrupted neuronal activity (shown in 
red and black) in the parameter plane of coupling strengths 
(cs, ea)-The maximum voltage from all neurons , averaged over 
10 random realizations of the network of 10 neurons is used 
as the index and specific nature of activity is verified from 
time series of a few representative cases. Black regions cor- 
respond to regions with total absence of activity and the red 
regions correspond to regions of sub-threshold activity. Yellow 
regions correspond to normal neuronal activity which can in- 
clude unsynchronized bursts and spikes, synchronized bursts, 
and synchronized spikes for different values of Ca, and Cs 







m 




i 


m 


w 




i 




i 



1 000 1 040 1 080 

Time 



FIG. 4: Revival of neuronal activity by preventing leakage 
of induced activity of A/3. A, Voltage of a neuron recovering 
from sub-threshold state with Cs = 0.5, Ca = 3.5 .B, Voltage 
of a neuron recovering from silent state with Cs = 1.0 and 
Ca = 2.0. In both cases, k is reduced to 1% of its initial value 
at the 1000*^ time unit. 
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